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The results of the muon-spin rotation experiments on BCS superconductors Nbl^+x (x = 0.2, 
0.34) are reported. Both samples, studied in the present work, exhibit rather broad transitions to 
the superconducting state, suggesting a distribution of the volume fractions with different transi- 
tion temperatures (T c )'s. By taking these distributions into account, the dependence of the inverse 
squared zero-temperature magnetic penetration depth (A^" 2 ) on T c was reconstructed for tempera- 
tures in the range 1.5 K< T c < 8.0 K. A^" 2 was found to obey the power law dependence Xq 2 oc T^' 1 
which appears to be common for some families of BCS superconductors as, e.g., Al doped MgB2 
and high-temperature cuprate superconductors as underdoped YBa2Cu307-<s. 
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I. INTRODUCTION 

The search for universal correlations between physi- 
cal variables, such as transition temperature, magnetic 
field penetration depth, electrical conductivity, energy 
gap, Fermi energy etc. may provide hints towards a 
unique classification of different superconductors. In par- 
ticular, establishing of such correlations may help to un- 
derstand the phenomenon of superconductivity, that is 
observed now in quite different systems, such as simple 
metals and alloys, fullerenes, molecular metals, cuprates, 
cobaltites, borides etc. Among others, there is a corre- 
lation between the transition temperature (T c ) and the 
inverse squared zero-temperature magnetic field pene- 
tration depth (Aq 2 ), that generally relates to the zero- 
temperature superfluid density (p s ) in terms of p s oc 
Aq" 2 . In various families of underdoped high-temperature 
cuprate superconductors (HTS)'s there is the empirical 
relation T c oc p s oc Aq 2 , first identified by Uemura et 
al**! It was recently shown, however, that for HTS's with 
highly reduced T c 's the direct proportionality between T c 
and Aq -2 is changed to a power law kind of dependence 
with {T c ) n oc Aq 2 (n is the power law exponent). In 
experiments on highly underdoped YBa2Cri307_5 Zuev 
et al£ obtained n = 2.3, Liang and coworkers reported 
n = 1.6, 4 while Sonier et al^ found n = 2.6 — 3.1. In 
molecular superconductors Pratt and Bhmdcll obtained 
n = 2/3j£ From the theoretical point of view, it was 
shown that in systems obeying 2D or 3D quantum su- 
perconductor to insulator transition, n = 1 or n = 2, 
respectively 

It should be emphasized, however, that the relation 
between T c and Aq~ 2 is not yet established for BCS 
superconductors and still awaits to be explored. A 
good candidate to search for such a relation could be 
NhE^+a;. The superconductivity in NhE^+a;, similar to 
MgB2, is most likely mediated by phonons. It is con- 
firmed in nuclear magnetic resonance (NMR)^ 1 - and tun- 



nelling experiments ; 12 ^ 3 as well as by recent calculations 
of the elastic properties^ Moreover, muon-spin rotation 
(/iSR) experiments suggest that the superconducting gap 
is isotropic^ As is shown in Refs. [l6| and H3 the su- 
perconductivity in NbB2 can be induced by either in- 
creasing boron or decreasing niobium content, while the 
parent NbB2 compound is not superconducting at least 
down to 2 K. The transition temperature was found to 
reach the maximum value of 9.2 K and 9.8 K for Nbo.gB2 
and NbB2.34, respectively] 16 i 17 This offers a possibility to 
study the relation between T c and Ag~ 2 as a function of 
boron and/or niobium content. In the present study the 
temperature dependence of the magnetic field penetra- 
tion depth was measured for two NhE^+z samples with 
x = 0, 2 and 0.34 by means of transverse-field muon-spin 
rotation technique. It was found that in both samples 
the distribution of the superconducting volume fractions 
with different T c 's can be well approximated by a Gaus- 
sian distribution. The mean values of the superconduct- 
ing transition temperature (T™) and the width of the 
distribution (AT C ) were found to be T c m = 6.02(3) K, 
AT C = 0.96(2) K for NbB 2 . 34 at p H = 0.1 T, and 
T™ = 3.40(4) K, AT C = 1.06(2) K for NbB 2 . 2 at 
PqH = 0.05 T. Within the model, developed for a gran- 
ular superconductor of moderate quality, we reconstruct 
the dependence of the zero-temperature superfluid den- 
sity p s oc Aq" 2 on the transition temperature T c . It was 
found that in the range of 1.5 K< T c < 8.0 K A^ 2 follows 
a power law dependence, rather than a linear dependence 



reported by Takagiwa et al.^ with A oc T, 



^3.1(1) 



The 



value of the power law exponent 3.1(1) agrees rather well 
with n = 2.6 — 3.1 reported by Sonier et al^ for under- 
doped HTS's YBa 2 Cu307-,5. 



The paper is organized as follows. In Sec. Ill Al we de- 
scribe the model used to obtain the temperature depen- 
dence of the magnetic field penetration depth for a gran- 
ular superconductor having a certain distribution of the 
superconducting volume fractions. The distributions of 
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the local magnetic fields, calculated within the framework 
of this model, are presented in Sec. Ill Bl In Sec. IIIII we 
describe the sample preparation procedure and details of 
the muon-spin rotation and magnetization experiments. 
Sec. IIVI comprises studies of the magnetic penetration 
depth in NDB2.2 and NDB2.34 superconductors. The con- 
clusions follow in Section fVl 



II. THEORETICAL BACKGROUND 

A. Magnetic penetration depth in a granular 
superconductor of moderate quality 

In this section we describe the model applied to calcu- 
late temperature dependence of the magnetic field pene- 
tration depth A in a granular superconductor of moderate 
quality by using the /iSR data. We based on a general 
assumption that A can be obtained from the second mo- 
ment of the local magnetic field distribution [-P(-B)] inside 
the superconducting sample in the mixed state measured 
directly in /iSR experiments . The model uses following 
assumptions: (i) The superconducting grains are decou- 
pled from each other, (ii) Each i— th grain is a super- 
conductor with a certain value of the transition temper- 
ature (T*) and the zero-temperature magnetic penetra- 
tion depth (Ag). (iii) The zero-temperature supercon- 
ducting gap (Ao) scales with T c in agreement with the 
well-known relation Aq/ ksT c = const^ implying that 
the ratio R = A l /ksT* is the same for all the grains. 
Note that a linear decrease of both superconducting en- 
ergy gaps (A CT and A 71 ') with decreasing T c was observed 
recently by Gonelli et aL— in Mn doped MgB2 . The sim- 
ilar linear Ao vs. T c scaling was reported by Khasanov 
et for RbC^Og BCS superconductor. 

Let us first define variables and functions used within 
the model. Function u(t) describes the distribution of 
the superconducting volume fractions with different tran- 
sition temperatures T c 's. It is defined so that the vol- 
ume fraction of the sample, having transition tempera- 
tures in the range between T % c and T-!, is obtained as 

J T i oj(t)dt. The distribution of the local magnetic fields 

in the i— th grain is described by the function P l (B), 
which for ideal superconductor has rather asymmetric 
shape [see Fig. Q] (a)]. Function /(£) describes depen- 
dence of the inverse squared magnetic penetration depth 
A" 2 at T = on T c [(Aj,)~ 2 = /(T*)]. The dependence 
of A* on temperature was assumed to follow the standard 
equation for weak coupled BCS superconductor [see, e.g., 
Eq. (2-111) in Rcf.[l|: 

[A*(T)r 2 = (A )" 2 S (T, A*) = /(T C >(T, k B R ■ T*), (1) 

with the temperature dependent part 

s(T, A n ) = 1 + 2/ ^— — dF, 

K °' J A i (r) \dEj ^/E 2 — A i (T) 2 



Here, F = [1 + cxp(_E/fcsT)]~ 1 is the Fermi function, 
A Z (T) = A A(T/T^) represents the temperature depen- 
dence of the energy gap, and A = 2RksTl denotes 
the zero temperature value of the superconducting gap. 
For the normalized gap A(T/T C ) the values tabulated in 
Ref. Hi] were used. Finally, the temperature dependence 
of the total second moment of the local magnetic field 
distribution in the superconducting sample can obtained 
as: 

2 rco roc 

(AB 2 ) tot = _ = / / cu(t)P{B',t)(B'-B m ) 2 dtdB'. 
% Jt Jo 

(2) 

Here, 7^ = 2ir x 135.5342 MHz/T is the muon gyromag- 
netic ratio, a 2 is the second moment of the ^SR line and 
B m is the mean internal field inside the superconducting 
sample. 

As will be shown later, in transverse-field /iSR ex- 
periments one obtains directly: (i) the distribution of 
the superconducting volume fractions with different T c 's 
[a->(T c )], (ii) the temperature dependence of the mean in- 
ternal field B m (T), and (iii) the temperature dependence 
of the second moment of the ^SR line a 2 (T) . By substi- 
tuting w(T c ) and B m {T) to the Eq. © and then fitting 
it to the experimental cr(T) data, one should be able to 
obtain the distribution of the zero-temperature super- 
fluid density p s oc Ag~ 2 as a function of the transition 
temperature T c and the ratio R = Ao/fcsT c . In order 
to do that one needs, however, to calculate P l (B) dis- 
tributions which depend on the applied field (B ex ), A 4 , 
and the coherence length (£') in a nontrivial way. This 
makes fit of Eq. to the experimental data extremely 
difficult. The analysis can be simplified by assuming that 
each P l (B) follows the Gaussian distribution and is de- 
termined, therefore, by the only two parameters: the sec- 
ond moment of the Gaussian line 

(AB 2 ) = (a/ 1 ,) 2 = G 2 (b)-(\)-\ 

and the internal field of the i— th grain (B l ). Here, 
b = B/B C 2 is the reduced magnetic field (B C 2 is the 
second critical field) and the G(b) is the proportional- 
ity coefficient between 17/7^ = \J (AB 2 ) and A -2 which 
can be obtained by means of Eq. (13) from Ref. as: 

G(b) = 0.172-^(1 - 6)[1 + 1.21(1 - Vb)% (3) 

($0 is the magnetic flux quantum). Here we also take 
into account that within Ginzburg-Landau theory £ = 
\J &o/2irB C 2- The internal field of the i— th grain was 
calculated within the London approximation modified by 
Brandt:— 

~ B ex - (1 - D m )^f(T*)a(T, k B RTl) Hg{b m W) 
where 

g(b) = 1 + ^-^(0.357 + 2.8906 - 1.581& 2 ), 
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V = B l /B l c2 , and b m = B m /B 



<:2 



B l c2 and D l are the sec- 



ond critical field and the demagnetization factor of the 
i — th grain, and D m and B^ are the mean values of the 
demagnetization factor and the second critical field, re- 
spectively. In the last part of Eq. ((4]) , we also replaced D % 
and b l by their average values D m and b m (it is clear, that 
in the intermediate range of fields and for temperatures 
not very close to T c , b % ~ b m ). 

The advantage to use Gaussian P % (B) is that the con- 
tribution of the i— th line to the total second moment can 
be obtained analytically as : 23 > 24 



(AB y = (a*/ lfi y + (B* 



B 



rn\2 



G 2 {b m )[\ l {T)]- i + (B l - B m f. 
By substituting it in Eq. (J5J one gets: 



(AB 2 ) tot = — 



oj{t){AB 2 )dt. 



(5) 



(6) 



Later on we are going to use this equation in order to 
fit the experimental /xSR data. We want to remind that 
the total second moment (AB 2 ) tot obtained by means 
of Eq. ijHJ) is determined by: (i) Function f(T c ) that 
describes distribution of Aq 2 as a function of T c [see 
Eq. (JTJ>] . (ii) Function u>(T c ) describing distribution of 
the superconducting fractions with different T c 'a. (iii) 
Gap to T c ratio R = 2A Q /k B T c [sec Eq. ©]. (iv) The 
mean values of the demagnetization factor D m and the 
reduced magnetic field b m [see Eq. (Q| and ©]. 

Note that in real experiments on polycrystallinc 
samples with sharp transition to the superconducting 
state, P(B) often obeys Gaussian distribution (see, e.g., 
Refs. [23| and l25h. The reason is that the vortex lattice 
within the small superconducting grain does not remain 
regular. It is distorted due to effects of pinning inside the 
grain, as well as, near the edges of the grain due to close- 
ness of vortexes to the surface. As is shown by Brandt^ 
even small pinning leads to substantial smearing of the 
characteristic features of the ideal P(B) line and, in a 
case of intermediate pinning, to the Gaussian shape of 
P(B). We should emphasize, however, that even in a 
case of distorted vortex lattice, the second moment of 
P(B) is still a good measure of A<2& 



B. Magnetic field distribution in a granular 
superconductor 

In this section we are going to simulate the internal 
magnetic field distribution P{B) within the above de- 
scribed model. As a first step it is necessary to choose a 
theoretical model describing the spatial variation of the 
local internal magnetic fields in the vortex lattice [-B(r)], 
from which the distribution of the fields in the i— th grain 
can be obtained as: 



P\B) 



J5{B~B') dA(B') 



(7) 



Here, dA(B') is an elemental piece of the vortex lat- 
tice unit cell where the magnetic field is equal to B' 
and the unit cell has a total area of / dA(B'). B(r) 
was calculated by using an iterative method for solving 
the Ginzburg- Landau equations developed by Brandt *22 
This method allows to accurately determine B(r) for ar- 
bitrary^,^ = A/f and the vortex lattice symmetry (see 
Refs. \2% [271. and [23 for details). It was also assumed that 
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FIG. 1: (Color online) P(B) distributions calculated by 
means of Eq. JTJJJ at T c m = 6 K, Ao(T c m ) = 60 nm, k = 1.67, 
and B ex = 0.1 T for the following values of AT C : 0.0 K 
(a) 0.2 K (b), 1.0 K (c), and 2 K (d). The dashed green, 
solid blue, and doted red lines in (b), (c), and (d) represent 
the coiT^P^B) term for T c = T" 1 - AT C , Tl = T c m , and 
T l c = T c m + AT C , respectively. See text for details. 

the distribution of the superconducting volume fractions 
with different T c 's is described by a Gaussian distribu- 
tion: 



1 



AT C V2^ 



exp 



2AT 2 



(8) 



(T™ and AT C are the mean value and the width of the 
distribution, respectively), A^ 2 follows the power law: 



(9) 



with the exponent n = 2, and k % = A J /C' = const^- As is 
mentioned already in the introduction, the power law de- 
pendence of A^~ 2 on T c is observed for various HTSii 2 -^^ 
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and molecular superconductors^ as well as obtained the- 
oretically for materials having 2D or 3D quantum super- 
conductor to insulator transition.— iSuLiii 

As initial parameters for calculations we took T" 1 = 
6 K, B ex = 0.1 T, D m = 1/3, A (T c m ) = 60 nm and « = 
1.67. The resulting field distribution P(B) was obtained 
as: 



JV 



(10) 



Calculations were done for N = 60 in the region ±3AT C 
around T™. Figure [1] shows P(B) distributions for 
AT C = 0.0 K, 0.2 K, 1.0 K and 2.0 K obtained by means 
of Eq. (|T0|) . The lines in (b), (c), and (d) represent the 
u>{Tl)P % {B) term for T* = T™ — AT C (dashed green line), 
T l c = T™ (solid blue line), and T* = T™ + AT C (dotted 
red line). It is seen that the shape of P(B) changes dra- 
matically with increasing width of the lu(T c ) distribution. 
For small AT C (when transition to the superconducting 
state is very sharp) P{B) is asymmetric with the highest 
weight around the point corresponding to the so called 
"saddle point" field B sa d [see Figs. [1] (a) and (b)]. It is 
also seen that all the characteristic features of P(B) at 
minimum (B rn i„), maximum (B max ), and "saddle point" 
fields are smeared out [see Fig.Q](b)]. Note that the sim- 
ulated P(B) presented in Fig. [1] (b) looks very similar 
to what is observed in fj.SK experiments on high-quality 
single crystals ! 30 ' 31 With a further increase of AT C the 
P(B) distribution becomes rather symmetric and, finally, 
asymmetric again, but now with the maximum weight 
around the external field B ex = 0.1 T and a very long 
tail at lower fields. It is interesting to note, that the 
parts of the sample with the lowest T c 's are responsible 
for the peak appearing slightly below the external field, 
as is seen in Fig.Q](c) and (d). The reason is the decrease 
of P' l (B) width and the shift of B l towards to B ex with 
increasing A 4 . 

To summarize, in Sec. [TT] we described the model al- 
lowing to obtain the distribution of the internal magnetic 
fields in a granular superconducting sample of moderate 
quality and calculated the second moment of this dis- 
tribution. Within the framework of this model we also 
simulated P(B) profiles for materials with different width 
of the superconducting transition AT C . It is remarkable, 
that already small AT C leads to smearing of the char- 
acteristic features of P(B) distribution near B m i n , B sa d, 
and B max characteristic fields. Even though this result is 
quite predictable, it has an important impact, since pre- 
viously smearing of P{B) was ascribed entirely for the 
effects of pinning (see, e.g., Ref. l3lh. 



III. EXPERIMENTAL DETAILS 

Details of the sample preparation for NhB^+z can be 
found elsewhere^. Both NhB2.2 and NbB2.34 samples, 
studied in the present work, were fine powders with the 
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FIG. 2: (Color online) Temperature dependences of the field- 
cooled magnetization (Mfc) for NDB2.2 and NbB2.34 samples 
measured at /j.oH = 0.5 mT. The solid lines are the best linear 
fits to the steepest part of Mfc{T) curves and M = 0. 



average grain size of the order of few microns. The field- 
cooled 0.5 mT magnetization (Mpc) measurements for 
NbB2.2 and NbB2.34 samples were performed by using 
a SQUID magnetometer. The corresponding M F c(T) 
curves are shown in Fig. [21 It is seen that the super- 
conducting transitions are rather broad indicating that 
both samples are not particularly uniform. This also im- 
plies that the superconducting critical temperatures may 
be evaluated only approximately The middle-points of 
transitions correspond to ~ 9.16 K and ~ 3.60 K, while 
linear extrapolations of Mp q (T) curves in the vicinity of 
T c to M = result in 9.76 K and 4.79 K for NbB 2 . 2 and 
NbB 2 .34, respectively (see Fig. [5]). 
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FIG. 3: (Color online) The muon-time spectra (a), differ- 
ence between the two-Gaussian fit and experimental data 
(b), and internal field distributions (c) for NbB2.34 sample at 
T = 1.7 K after field cooling in magnetic field of 0.1 T. The 
lines in (a) and (c) represent the best fit with the Gaussian 
line-shapes. See text for details. 

The transverse-field /iSR experiments were performed 
at the 7rM3 beam line at the Paul Scherrer Institute (Vil- 
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ligen, Switzerland). The powders of NDB2.2 and NDB2.34 
were cold pressed into pellets (12 mm diameter and 2 mm 
thick). The samples were field cooled in fi H = 0.05 T 
(NDB2.2) and 0.1 T (NI3B2.34), applied perpendicular to 
the flat surface of the pellet, from above T c down to 1.7 K. 
The fields 0.05 T and 0.1 T were chosen in order to per- 
form measurements at the same reduced magnetic field 
B/B c2 ~ 0.4 (B c2 is the upper critical field) M The /iSR 
signal was observed in the usual time-differential way by 
counting positrons from decaying muons as a function of 
time in positron telescopes. The time dependence of the 
positron rate is given by the expression^ 

=N -e- t ^[l + AP(t)] + bg, (11) 

where Nq is the normalization constant, bg denotes the 
time-independent background, r M = 2.19703(4) x 10 -6 s 
is the muon lifetime, A is the maximum decay asymmetry 
for the particular detector telescope (A ~ 0.18 — 0.19 
in our case), and P(t) is the polarization of the muon 
ensemble: 

P(t) = J P(B) cos(7 M Bt + (j))dB . (12) 

Here, <£> is the angle between the initial muon polariza- 
tion and the effective symmetry axis of a positron de- 
tector. P(t) can be linked to the internal field distribu- 
tion P(B) by using the algorithm of Fourier transform^ 
The P(t) and P(B) distributions inside the NbB2.34 at 
T = 1.7 K after field-cooling in a magnetic field of 0.1 T 
are shown in Fig. [3] (a) and (c). The P(B) distribu- 
tions were obtained from measured P(t) by using the 
fast Fourier transform procedure based on the maximum 
entropy algorithm^ In order to account for the field dis- 
tribution seen in Fig. [3] (c), the /iSR time-spectra were 
fitted by two Gaussian lines: 

2 

P(t) =Y, A i exp(-cr?t 2 /2) cos{^Bit + cp) . (13) 

i=l 

where A{, ai, and Bi are the asymmetry, the Gaus- 
sian relaxation rate and the first moment of the i-th 
line. At T > 8 K for NbB 2 . 34 and T > 4.5 K for 
NbB 2 .2 the analysis is simplified to the single line only 
with a nm ~ 0.3 fis^ 1 resulting from the nuclear mo- 
ments of the sample. Eq. (fl~3|) is equivalent to the field 
distribution)^ 

P ^ B ) = Tm E - -P [ ^2 j ■ (14) 

The solid line in Fig. [3] (a) represent the best fit with 
the two-Gaussian lines to the muon-time spectra. The 
corresponding P(B) is shown in Fig. [3] (c). It should be 
mentioned that the two-Gaussian fit can satisfactory de- 
scribe the experimental data. For both samples and in 
the whole range of temperatures the normalized Xnorm' s 
were found to be close to the unity implying a good qual- 
ity of fits. 



IV. RESULTS AND DISCUSSIONS 

Results for NbB2.2 and NbB2.34, obtained from the fit 
of experimental data by means of Eq. Q13p , are summa- 
rized in Fig. 0] Distributions of the local fields, at tem- 
peratures marked by the arrows, are shown in the upper 
panels of Figs. 0] (a) and (d). The broad lines reflect 
contributions of the superconducting parts of samples. 
Indeed, their asymmetries A sc [Figs. 3] (a) and (d)] and 
relaxation rates a sc [Figs. [4] (c) and (f)] increase while 
first moments B sc [Figs, (b) and (e)] decrease with de- 
creasing temperature, as expected for a superconductor 
in a mixed state (see, e.g., Ref. [35]). The narrow lines 
describe contributions from parts of the samples being 
in the normal state. The slight shift of these lines to 
higher fields [Figs. Q] (b) and (e)] and the small increase 
of relaxations [Figs. SI (c) and (f)] are associated with 
the diamagnetism of the superconducting grains leading 
to increase of the local fields in the nonsuperconduct- 
ing parts of the sample. Such behavior is often observed 
in samples with less than 100% superconducting volume 
fraction (see, e.g., Refs. l2~il and [36h . 

In Fig. [5] we plot temperature dependences of the su- 
perconducting components A sc and cr sc . The a sc values 
presented in Fig. [5] are corrected to the nuclear moment 
contribution a nm which was subtracted in quadrature 
(see e.g. Ref.l24h. It is seen [Fig.[S](a)] that superconduc- 
tivity does not disappear abruptly. The superconduct- 
ing fraction decreases continuously from their maximum 
value to zero with temperature rising from 4 K to 8 K for 
NbB 2 . 34 and from 1.5 K to 5.5 K for NbB 2 . 2 . The solid 
lines in Fig. \5\ (a) correspond to fits of A SC (T), assum- 
ing that the distribution of the superconducting volume 
fractions with different T c 's [w(T c )] follows the Gaussian 
distribution [sec Eq. l[8"])] and, consequently, 

A SC (T) = A™ f°° uj(t)dt . (15) 
Jt 

Here, A™ is the mean value of the superconducting asym- 
metry at low temperatures. The fits yield T™=6.02(3) K, 
AT C =0.96(2) K, and A™=0.170(5) for NbB 2 . 34 , and 
T c m =3.40(4) K, AT C =1.06(2) K, and A™=0.128(7) for 
NbB2.2- By taking into account that the total asymme- 
tries (A sc + An) were found to be ~ 0.19 for NbB 2 .34 and 
~ 0.18 for NbB 2 . 2 [dotted lines in Figs.SJa) and (d)], the 
superconducting volume fractions were estimated to be 
-85 % and -70 % for NbB 2 . 3 4 and NbB 2 . 2 , respectively. 

The samples used in our /iSR experiments were cold 
pressed powders. In this case the individual grains are 
expected to be only weakly coupled. This implies that 
the model, described in Sec. Ill Al can be applied for the 
particular NbB2.2 and NbB2.34 samples studied in the 
present work. All the variables and functions needed for 
such calculations as, e.g., the distribution of the super- 
conducting volume fractions with different T c 's [uj(T c ) - 
dashed red and blue lines in Fig. [5] (a)], the temperature 
dependences of the square root of the second moment 
Kc(T) = j^(ABi)(T), sec Figs. @ (c) and (f)] and 
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FIG. 4: (Color online) The temperature dependences of asymmetries (A ac , A n ) (a)/(d), internal fields (B ac , B„) (b)/(e), 
and the muon-spin depolarization rates (<r sc , tj„) (c)/(f), obtained from the fit of jj.SK time spectra by means of Eq. (|13[1 for 
NbB2.34/NbB2.2 after field-cooling in a magnetic field of 0.1 T/0.05 T. The indexes "sc" and "n" denote the normal and the 
superconducting state components. The dotted lines in (a)/(d) and (b)/(e) represent the total muon asymmetry (A sc + A n ) 
and the external field value (B ex ), respectively. In the upper panels of (a) and (d) the P(B) distributions at temperatures 
marked by the arrows are shown. The red and the blue dotted lines represent the superconducting (broad) and the normal 
state (narrow) components obtained from the fit by means of Eq. (|13[) . The solid green line corresponds to the sum of these 
two components. 



the mean field [B SC (T) = B m (T), see Figs. H (b) and 
(e)] are directly obtained in the experiments. The only 
unknown parameter is the demagnetization factor D m 
which enters Eq. We should mention, however, that 
it is not correct to assume that D m is equal to the some 
average demagnetization factor for all the grains as, e.g., 
D m = 1/3 for grains of spherical symmetry. The reason 
for that is the following. The individual grains have dif- 
ferent internal magnetic fields and, therefore, will show a 
different diamagnctism. As a consequence, the internal 
field in the particular grain is determined by the diamag- 
nctism of the grain itself, by the demagnetization field 
of the whole sample and by the fields from the other su- 
perconducting grains surrounding it. This problem was 
already discussed by Weber et dl. in Ref. [23|. According 



to their calculations for the grains of spherical symmetry 
the factor 1 — D m in Eq. Q should be replaced with [see 
Eq. (36) in Ref.Hl]: 

1 - D m ~ 2/3 - (DP - 1/^Vp/VG- (16) 

Here, D p is the demagnetization factor of the whole sam- 
ple, and rj p and rjc are the effective volume density and 
the x-ray density of the sample, respectively. For the 
experimental geometry (thin disk in perpendicular mag- 
netic field) D p ~ 1. Assuming now that the density of 
the pellet is twice as small as the x-ray density of the 
material (rjcMp = 2) we get 1 — D rn ~ 0.3. At this 
stage we are not going to estimate the value of 1 — D" 1 
more precisely. As is shown below, it can be obtained self 
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FIG. 5: (Color online) (a) Temperature dependences of the 
superconducting asymmetries (^4 SC ) for NDB2.34 and NbB2.2- 
The solid lines are the fits of A ac (T) data by means of Eq. ()15|l . 
The dashed lines represents the Gaussian distributions of the 
superconducting volume fractions with different T c 's [cj(T c )]. 
(b) Temperature dependences of a sc corrected to the nuclear 
moment contribution a nm for NbB2.34 and NbB2.2- The open 
and filled small black circles are the fit of er sc (T) by means of 
Eq. (j6j). See text for details. 



consistently from the fit of Eq. ([6]) to the experimental 
data. 

The demagnetization factor D m enters Eq. ^ via the 
term (AS 2 )* that, in its turn, is a sum of G 2 (6) [A* (T)]" 4 
and (B l - B m f [see Eq. ©]. The term (B l - B m ) 2 , 
which depends on D m , is responsible for the correction 
to (AB 2 ) 1 appearing due to the shift of the internal field 
of the i— th grain from the mean internal field B m [see 
Eq. (J5J)]. It is clear, that during fit of Eq. to the 
experimental data, contribution of (B l — B m ) 2 term to 
(AB 2 ) tot needs to be minimized (B m should become the 
first moment of the resulting P{B) distribution). We 
used, therefore, an iterative approach. During the fit, 
both a sc (T) data sets for NbB2.2 and NbB2.34 samples 
were fitted simultaneously. The dependence of Ao on the 
transition temperature was assumed to be described b 
the power law A^ 2 = (K ■ T c ) n [see Eq. © and Refs " 
0@@QIlS[i3] ■ I n a nrs t step, 1 — D m was assumed to 
be equal to 0.3 (see above) and Eq. ([6]) was fitted to 
the data with the proportionality factor K, the power 
law exponent n, and the gap to T c ratio R = Aq/UbTc 
as free parameters. In a second step, the values of K, n, 
and R, obtained in the step one, were substituted back to 
Eq. © and the second term (B l — B m ) 2 entering (AB 2 ) 
in Eq. ^ was minimized with the only free parameter 
D m . Then the whole cycle was repeated by using as 



initial parameter the newly obtained D rn value. After 3 
iterations the fit already converges. The fit yields K = 
1.24(3) • 1CT 2 nm-^K" 1 , n = 3.1(1), R = 1.68(3) and 
1 - D m = 0.24(2). The open and the filled black circles 
in Fig. [5] (b) represent the result of the fit of Eq. [6] to the 
data. 

Three important points emerge: 

(i) The value of R = A /k B T c = 1.68(3) obtained from 
the fit is very close to the weak-coupling BCS value 1.76 
and Ao/ksTc = 1.55 obtained by Kotegawa et al^ in 
n B NMR experiments. It is, however, much smaller than 
A /k B T c = 2.15-2.25 reported by Ekino et alM in tun- 
neling experiments. 

(ii) The value of 1-D m = 0.24(2) was found to be rather 
close to 0.3 roughly estimated from Eq. (fT6f . 

(iii) For each particular data set the Aq -2 vs. T c depen- 
dence can be reconstructed only for T c 's in the range 
T[ < T c < Tj} (T' c l and T\ denote the temperatures at 
which the superconducting fraction achieves the maxi- 
mum value and vanishes, respectively). Figure O (a) re- 
veals that the corresponding regions are from ~ 4 K to 
- 8 K and from ~ 1.5 K to ~ 5.5 K for NbB 2 . 34 and 
NbB2.2, respectively. Bearing in mind that fit by means 
of Eq. @ was performed for both NbB 2 .34 and NbB 2 .2 
data sets simultaneously, we can thus conclude that the 
relation A^~ 2 oc Tc is valid at least for temperatures 
in the region from ~ 1.5 K to ~ 8 K. 

In Fig. [S] we plot T c vs. \^ 2 dependence obtained 
from the fit of NbB2.34 and NbB2.2 data. The grey 
area represents the region where A^ 2 = [(1.24 ± 0.03) • 
10" 2 • T c } 3 • 1±0 - 1 . The dotted line corresponds to Aq 2 = 
[1.24 ■ 10~ 2 • Tc] 3 - 1 . The dashed blue and the dot-dashed 
red lines represent the range of the transition tempera- 
tures where T c vs. Aq 2 dependences were reconstructed 
for NbB 2 .2 and NbB 2 .34, respectively. For clarity, these 
lines are shifted by 0.05 K above and below the central 
(1.24 • 10~ 2 • Tc) 3 - 1 line. We also include in this graph 
data points for NbB 2 + K (x — 0.0, 0.01, and 0.1) from 
Ref . Il5l . It is seen that these samples have approximately 
twice as high transition temperatures as one would ex- 
pect from the obtained T c vs. Aq 2 dependence. On the 
base of our results, we can argue that samples studied 
in Ref. [TH also have distributions of the superconduct- 
ing volume fractions, similar to what is observed in the 
present study. Without taking into account this distribu- 
tions, fit of <J SC (T) data can lead to a substantial overes- 
timate of the superconducting transition temperature T c . 
As an example, the star in Fig. [6] represents result of the 
fit of Eq. ([TJ) to cr sc (T) NbB 2 .2 data when all measured 
points arc taken equally into account. 

Now we are going to comment shortly the observed T c 
vs. Aq 2 dependence. Recently it was shown that the fa- 
mous "Uemura" relation, establishing the linear propor- 
tionality between T c and Aq 2 ^' 2 - does not hold for highly 
underdoped HTS's^ Indeed, for YBa 2 Cu 3 7 -5 Zuev 
et al£ observed Aq 2 cx T 2 ' 3 ^' and found that this power 
law is in fairly good agreement with YBa2Cu307_,5 data 
in the whole doping range. For the similar compounds 
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FIG. 6: (Color online) Dependence of the transition temper- 
ature (T c ) on the inverse squared zero-temperature magnetic 
penetration depth (Xq 2 ) for NbB 2 +;c. The grey area repre- 
sents the region where Aq 2 = [(1.24 ± 0.03) • 10~ 2 • T c ] 3 1±01 . 
The dotted line corresponds to A c 7 2 = (1.24- 10" 2 -T c ) 31 . The 
dashed blue and the dash-dotted red lines represent the range 
of transition temperatures where T c vs. Xq 2 dependences were 
reconstructed for NbB 2 . 2 (1.5 K< T c < 5.5 K) and NbB 2 . 3 4 
(4 K< T c < 8 K). The lines are shifted by 0.05 K above and 
below the central (1.24 ■ 10 -2 ■ Tc) 3 ' 1 line for clarity. The star 
is the T c vs. Xq 2 point obtained from the fit of a ac (T) data 
for NbB 2 . 2 by means of Eq. (JTJ when all the measured points 
are equally included in the fit. The open circles are the data 
points from Ref. [l5|. Values of A(0) -2 were obtained from 
(Tsc(0) measured in a field 0.1 T and # c2 (0) (see Ref. (3) by 
using Eq. (13) from Ref.[H 



Sonicr et aL— obtained n = 2.6 — 3.1. Those values are 
very close to 3.1(1) obtained in the present study for 
NbB2+i; superconductor. The observation of the power 
law type of relation between the transition tempera- 
ture T c and the superfluid density in BCS supercondutor 
NbB 2+:E and their good correspondence to what is ob- 
served in HTS's points to close similarity between these 
materials. It, probably, comes from the fact that the in- 
crease of T c in NbB 2+:E , similar to HTS, is determined by 
increasing the charge carrier concentration. 15 ' 37 Indeed, 
the x-ray photoelectron spectroscopy experiments of Es- 
camilla and Huerta 3 - 7 - show that the increase of boron 
content leads to a decrease in the contribution of the Nb 
Ad states and increase in the contribution of the B 2p n 
states to the density of states at the Fermi level. 

It is important to emphasize that observation of cor- 
relation between T c and Aq 2 in BCS superconductors is 
not restricted to the particular NbB 2+:E system studied 
here. As an example, in Fig. [7] we plot T c /T™ ax as a 



function \q /{X 7 £ ax )- 2 for Al doped MgB 2 from Ref. 123. 
Here T™ ax is the maximum T c of a certain superconduct- 
ing family and X™ ax is the corresponding zero tempera- 
ture penetration depth. We also include on this graph 
points for YBa 2 Cu307_,5 from Ref. H|. It is seen that 
all these superconductors represent the very similar rela- 




♦ Mg„ Al B„ Ref. 38 

a 1-X X 2' 

★ YBa 2 Cu 3 7 5 , Ref. 5 



FIG. 7: (Color online) Dependence of the transition tem- 
perature normalized to the maximum transition tempera- 
ture of the superconducting family (T c /T™ ax ) on the normal- 
ized inverse squared zero-temperature magnetic penetration 
depth [\~ 2 /(X™ ax )~ 2 ]. Solid line corresponds to A^ 2 oc T c 3 1 
as obtained in the present study. Blue diamonds are the 
Mgi-a; Ala;B 2 data from Ref. l3a. Red stars are the data points 
for YBa 2 Cu307_i from Ref. [3. 



tions. The scaling relation for the BCS superconductors, 
reported here, and their agreement with what was ob- 
served in high-temperature cuprate superconductors 3 ^^ 
strongly suggests that there are some features of their 
electronic properties that are common, despite these ma- 
terials have quite different dimensionality, Fermi surface 
topology, symmetry of the superconducting order param- 
eter etc. 



V. CONCLUSIONS 

Muon-spin rotation studies were performed on BCS su- 
perconductors NbB2+a; (x = 0.2, 0.34). As a real space 
microscopic probe, fj.SK allows to distinguish between the 
superconducting and nonsuperconducting parts of the 
samples and determine the distributions of the supercon- 
ducting volume fractions with different T c 's. By using the 
model, developed for a granular superconductor of mod- 
erate quality, the dependence of the zero-temperature su- 
perfluid density p s oc A^ 2 on the transition temperature 
T c was reconstructed in a broad range of temperatures 
(1.5 K< T c < 8.0 K) revealing p s oc A„ 2 oc T c 31(1) . This 
dependence appears to be common at least for some fam- 
ilies of BCS superconductors as, e.g., Al doped MgB 2 
and high-temperature cuprate superconductors as, e.g. 
YBa 2 Cu 3 7 -5. 
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